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ABSTRACT 


L annular canister length, cm 


The Solar Dynamic Power Module being deve- 
loped for Space Station Freedom uses a eutectic 
mixture of LiF-CaF 2 phase-change salt contained 
in toroidal canisters for thermal energy storage. 
This paper presents results from heat transfer 
analyses of the phase-change salt containment can- 
ister. A two-dimensional, axisymmetric finite- 
difference computer program which models the 
canister walls, salt, void, and heat engine work- 
ing fluid coolant was developed. Analyses 
included effects of conduction in canister walls 
and solid salt, conduction and free convection in 
liquid salt, conduction and radiation across saLt 
vapor filled void regions and forced convection in 
the heat engine working fluid. Void shape, loca- 
tion, growth or shrinkage (due to density differ- 
ence between the solid and liquid salt phases) 
were prescribed based on engineering judgement. 

The salt phase change process was modeled using 
the enthalpy method. Discussion of results 
focuses on the role of f ree-convect ion in the liq- 
uid salt on canister heat transfer performance. 

This role is shown to be important for interpret- 
ing the relationship between ground-based canister 
performance (in 1-g) and expected on-orbit perform- 
ance (in micro-g). Attention is also focused on 
the influence of void heat transfer on canister 
wall temperature distributions. The large thermal 
resistance of void regions is shown to accentuate 
canister hot spots and temperature gradients. 
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NOMENCLATURE 


P 


density, g/cm 3 


A 


void surface element area, cm 2 or 
constant 


Stefan-Bol tzmann constant, 
5.67051xl0" 12 W/cm 2 -K 4 


B constant 

c specific heat, J/g-K 

e specific enthalpy, j/g 

F view factor 

g gravitational acceleration, cm/sec 2 

H PCM heat of fusion, J/g 

k thermal conductivity, W/cm-K 


Subscripts : 

i inner radius 

j void surface element 

k void surface element 

L liquid PCM 

m PCM me 1 1 

o 


outer radius 


2 


s 

solid PCM 

V 

void 

WS 

side wall 

w • 

canister wall 


1. INTRODUCTION 

Electrical power for NASA's proposed Space 
Station Freedom will be generated by photovoltaic 
solar arrays initially and later augmented with 
Solar Dynamic Power Modules (SDPM’s). The SDPM 
employs a concentrator to collect and focus solar 
energy into a cylindrical cavity heat receiver 
where it is converted to thermal energy. A frac- 
tion of the thermal energy is transferred to a 
circulating working fluid to operate a heat engine 
and produce electrical power. The remaining ther- 
mal energy melts a eutectic composition LiF-CaF 2 
Phase Change Material (PCM) contained in multiple 
canisters located concentrically around working 
fluid tubes running the length of the heat receiver 
cavity shown conceptually in figure 1. A single 
PCM containment canister is shown in figure 2. 

The PCM stores and releases thermal energy by 
undergoing phase change at its critical tempera- 
ture of 1040 K. This permits continuous operation 
of the heat engine during the substantial eclipse 
periods (up to 36 min) of Freedoms low earth 
orbit. The service life requirement for the heat 
receiver is 30 years. 

In previous work by the authors [l], computer 
programs were developed to analyze PCM contain- 
ment canister thermal performance using the 
enthalpy method. First, approximate one-dimensional 
analyses were conducted to estimate the effects of 
PCM void vapor heat transfer and liquid PCM free 
convection on canister thermal performance. Sec- 
ond, two-dimensional (r,z), or 2D (r,z), ax i sym- 
metric canister analyses were conducted without 
the presence of a PCM vapor void or liquid PCM 
free convection. Results from these analyses 
indicated that the presence of a one-dimensional 
void increased peak containment canister wall tem- 
peratures as much as 200 K. The presence of free 
convection lowered wall temperatures as much as 
170 K and lowered PCM temperature gradients by a 
factor of ~3 times in addition to increasing PCM 
melting rate by ~13 percent. However, two- 
dimensional results indicated that a large portion 
(about 50 percent) of PCM containment canister 
radial heat transfer occurred in the container 
metallic walls and not in the PCM. Therefore, 
it was postulated that the effects of a void and 
free convection would be much less dramatic when 
two-dimensional analyses were conducted. 

To quantify the effects of a void and free 
convection, the 2D (r,z) canister computer program 
described in [l], was extended to include a void 
model and a free convection model. Figure 3 
illustrates the PCM containment canister sche- 
matic problem geometry. This paper presents 
results from two-dimensional PCM containment can- 
ister numerical analyses which include a PCM 
vapor void and liquid PCM free convection. 


The PCM vapor void model includes a prescribed 
void shape and location. Void heat transfer 
occurs via uncoupled vapor thermal conduction and 
internal void surface radiation. Liquid PCM 
free convection is modeled by use of an enhanced 
liquid PCM thermal conductivity which is based 
on the Nusselt number (Nu) . Results are discussed 
in a comparative manner, highlighting significant 
differences in PCM containment canister tempera- 
ture and phase distributions that arise from the 
inclusion of a void and free convection. 

2. ANALYSIS 

The governing equation for energy redistribu- 
tion in the PCM and canister walls is conserva- 
tion of energy. The conservation equation was 
formulated using ’’the enthalpy method’’ as follows: 

= div (k7T) (1) 


where e is specific enthalpy given in joules per 
gram [2]. Specific enthalpy is coupled to temper- 
ature through the following set of constitutive 
equations : 
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( 2 ) 

Here, T m is the PCM melting temperature, H m is 
the PCM heat of fusion, and eg, cl, and c w are 
the specific heat values for the solid PCM, liquid 
PCM, and canister wall material, respectively. 

A thorough discussion of the governing equa- 
tions, boundary condi t ions , material properties, and 
numerical method used in this analysis are described 
in [ 1 ] . The PCM vapor void model is discussed 
below. 

2.1 PCM Vapor Void Model 


Due to receiver design requirements and PCM 
contraction during solidification, a fraction of the 
containment canister total volume consists of PCM 
vapor void at all times during the orbital cycle. 

The void volume fraction (VVF) yaries between 8 per- 
cent, when all PCM is liquid at the melting point 
(T m ) , to 22 percent, when all PCM is solid at T m . 
The small volume changes associated with thermal 
expansion of the PCM and containment canister 
walls are ignored. The void geometry selected was 
a cylindrical annulus which easily conformed to the 
cylindrical finite-difference element grid network. 
The void was placed adjacent to the canister outer 
wall, a location that will generate conservatively 
high canister wall temperature distributions (fig- 
ure 3). Void growth or shrinkage occur uniformly 
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across the PCM-void interface defined as r v . As 
PCM liquifies or freezes, r v increases or 
decreases, respectively, about 0.1 cm which changes 
void volume. The associated PCM growth or 
shrinkage is accommodated numerically by the combi- 
nation of variable radial grid size and a variable 
PCM computational domain. This procedure, known as 
the "combined grid element technique," is described 
in [1]. However, at the time of this writing, the 
PCM growth-shrinkage computer algorithm was not 
thoroughly checked-out and thus, was not incorpo- 
rated into these analyses. Therefore, as an engi- 
neering approximation, a constant 15 percent VVF 
was assumed. 

It has been shown that both void vapor conduc- 
tion and radiation are important elements of void 
heat transfer [ 1 ] . Since the void vapor pressure is 
small (7xi0 “ 3 torr at 1040 K ) , vapor mass is negli- 
gible ( 10 - 8 g). Therefore, the void vapor tempera- 
ture distribution can be determined by the steady 
state heat diffusion equation: 

<»> 

Preliminary void conduction calculations show 
that axial temperature gradients are small and can 
be ignored for engineering calculations. This elim- 
inates the second term in equation (3). The solu- 
tion to equation (3) has the familiar form 

T(r) = A ln( r ) + B , (4) 

with the constant A = {T(r 0 ) - T ( r v ) } / 1 n ( r 0 /ry) and 
the constant B = T(r 0 ) - ln(r 0 ){T(r 0 ) - T(r v )}/ 
ln(r 0 /r v ). Equation (4) can be evaluated at each 
axial void grid element to determine the void vapor 
temperature distribution. 

Void vapor thermal conductivity, k v , was esti- 
mated as in [ 2 ] using the kinetic theory of gases. 
Based on this theory, k v can be numerically evalu- 
ated using the following equation: 

k v - 1 . 457x10”® * T 1/2 . (5) 

In this equation, k v has the units W/cm-K and 
the temperature, T (in degrees K), is determined by 
equation (4). 

Void radiation heat transfer is calculated 
based on the assumptions that, ( 1 ) all void sur- 
faces are diffuse and gray, ( 2 ) PCM surfaces are 
opaque to all wavelengths of radiation, and (3) void 
vapor is a nonparticipating medium. With these 
assumptions, the governing equation set for void 
radiation heat transfer is: 



( 6 ) 


where k indexes from 1 to NRS [3]. In equa- 
tion (6), k and j are void surface element num- 
bers that take on all integer values between 1 and 
NRS . NRS is the total number of radiating surfaces 
in the void enclosure. The term is the 

Kronecker delta function equal to 1.0 for k * j 


and equal to 0.0 for k * j. Given that each sur- 
face element temperature, T » , eraittance, e j , and 
element- to-element view factor, Fj^j , are known, the 
surface element net radiative heat loss matrix, 
Qj/Aj, can be determined. These heat loss terms are 
then added to the energy balance equations for the 
appropriate finite-difference grid elements in the 
canister outer wall, canister side walls, and outer- 
mos t PCM . 

An emittance value of 0.52 was selected for 
canister walls which are fabricated of Haynes alloy 
188 (HA 188). This value is based on experimental 
measurements from diffuse (grit blasted) HA 188 test 
coupons for the temperature range 1000 to 1100 K. 

An emittance value of 0.6 was selected for PCM 
surfaces. This value is an estimate based on emit- 
tance data for similar dielectric materials in the 
temperature range of interest. Element- to-element 
view factors, F^j, were determined using existing 
closed-form view factor solutions and considerable 
view factor algebra. View factors can be recalcu- 
lated as void size varies during the melt-freeze 
cycle. 

3. RESULTS 

In the absence of applicable exact or analyti- 
cal solutions, previous numerical solutions, and 
experimental data, numerical consistency checks were 
performed to assess the validity of numerical solu- 
tions. A numerical check of canister model energy 
balance and void surface element view factor summa- 
tion was carried out for each computer run. For ail 
cases, an energy balance was maintained within 
0.0015 percent and all surface element view factors 
summed to 1.0 within machine accuracy. 

3 . 1 Canister Wi th Void 

Figure 4 illustrates canister temperature con- 
tour and PCM phase maps at four times (24.28, 

48.56, 69.80, and 81.94 min) during a 91 min orbital 
melt-freeze cycle. The PCM melting portion of the 
cycle occurs from time = 0 to ~55 min while PCM 
freezing occurs for the remaining portion of the 
cycle from time = ~55 to ~91 min. Initially, all 
PCM is solid at time *= 0 min. During the PCM 
melting period, the large void thermal resistance 
forces a large percentage of canister heat transfer 
to occur via canister walls. This is illustrated in 
figure 4(a) which shows high temperature gradients 
in the void region and isotherm normals (i.e., the 
direction of heat flow) generally aligned parallel 
to canister walls. Thus, energy absorbed in the 
canister outer wall diffuses around the void, down 
the side walls, and then into the PCM. As a conse- 
quence of this heat flow pattern, PCM melting 
occurs axially inward from'both side walls. By time 
- 48.56 min when ~90 percent of the PCM is liquid 
(figure 4(b)), heat transfer axially along the can- 
ister inner wall initiates PCM melting radially 
outward until all the PCM is liquified at time 
= ~55 min. 

During PCM freezing, the heat of fusion 
energy liberated is transferred to the engine work- 
ing fluid that cools the canister inner wall and to 
the canister outer wall where radiative heat loss to 
the receiver cavity occurs. Because the void acts 
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as a thermal insulator, much of the heat loss from 
the liquid PCM to the canister outer wall occurs 
via conduction in canister side walls. As a conse- 
quence of this heat flow pattern, PCM freezing 
occurs along the canister inner wall and along the 
void surface and the maximum side wall temperature 
exists at about the radial midpoint (figure 4(c)). 
Near the end of the orbital cycle at time 
= 81.94 min, ~90 percent of the PCM is frozen and 
the last remaining liquid PCM exists adjacent to 
void (figure 4(d)). Side wall radial temperature 
profiles have been reestablished with temperature 
increasing in the positive radial direction. This 
is due to a small radiative heat input to the canis- 
ter outer wall from the receiver cavity during the 
last "18 rain of the orbital cycle. 

Figure 5 illustrates the corresponding canister 
temperature contour and PCM phase maps for a 91 min 
melt-freeze cycle which includes free convection in 
the liquid PCM. Figure 5(a) is identical to fig- 
ure 4(a) since the liquid PCM Rayleigh number (Ra) 
is below the critical Ra and liquid PCM heat 
transfer is still dominated by conduction. Thus, 
no convective heat transfer enhancement takes place 
during the early part of the orbital cycle. At 
about 30 min into the cycle, the critical Ra is 
exceeded and the Nu begins steadily increasing 
from 1.0 to a value of 4.5 at the end of the melt- 
ing period ( —55 min). Free convection in the 
liquid increases the rate of PCM melting and 
decreases canister temperature gradients as shown 
in figure 5(b) where Nu = 3.047. During the PCM 
freezing portion of the cycle, liquid convective 
effects quickly die out and the Nu falls back to 
1.0 by "60 min into the cycle. Thus, for the major- 
ity of the PCM freezing period, liquid PCM heat 
transfer is dominated by conduction. Hence, fig- 
ures 5(c) and (d) are nearly identical to fig- 
ures 4(c) and (d) with the exception that slightly 
more liquid PCM exists at the times of comparison 
for the case with free convection, i.e., the values 
of MFL (PCM mass fraction liquid) are slightly 
greater . 

Figure 6 illustrates void radial heat trans- 
fer, Qvoid’ as a function of time during the orbital 
melt-freeze cycle. Qvoid is comprised of vapor con- 
duction and surface-to-surface radiation components. 
By convention, these components are taken as posi- 
tive if the resulting heat transfer is radially 
inward. During PCM melting, void heat transfer 
via radiation is about 3 times greater than that by 
vapor conduction and both components are positive 
and remain fairly constant. During PCM freezing, 
radiation is about 2 times greater than vapor con- 
duction and both components remain negative until 
"85 min into the cycle when all the PCM has 
frozen. The jump in the curves at "72 min is asso- 
ciated with the outer wall radiative heat flux 
boundary condition going from negative to positive. 
The fact that both components of void heat transfer 
remain negative from "72 to ~85 min has interesting 
implications for the canister heat flow pattern. 

For this time period, relatively warm PCM trans- 
fers heat radially outward across the void to the 
canister outer wall where it is then transferred 
back down the canister side wall and into the engine 
working fluid coolant (figures 4(d) and 5(d)). 


3.2 Canister Performance Comparison 

In this section, results from three canister 
analytical cases are compared: (1) without a void 

[l], (2) with a void, and (3) with a void and liquid 
PCM free convection. Comparison of results from 
cases 1 and 2 is intended to show the impacts of a 
void on canister heat transfer performance. Compar- 
ison of results from cases 2 and 3 is intended to 
show likely canister performance differences during 
ground testing (in 1-g) and flight operation (in 
micro-g) . 

Figure 7 illustrates the maximum canister wall 
temperature, T(r,z) max , throughout the 91 min orbi- 
tal cycle for cases 1 to 3. T(r,z) max occurs at 
the axial midpoint of the outer wall and does not 
vary in position with time. The introduction of a 
void increases T ( r , z ) max about 20 K during the 
PCM melting period (0 to "55 min) and decreases 
T(r,z) max about 10 K during the PCM freezing 
period (~55 to ~91 min). These changes are associ- 
ated with the large void thermal resistance which 
increases canister outer wall temperature gradients 
to conduct heat around the void. In addition, the 
void increases outer wall axial temperature gradi- 
ents by "100 percent and increases side wall radial 
temperature gradients by ~80 percent during the PCM 
melting period. 

The inclusion of liquid PCM free convection 
reduces T(r,z) raax 0 to 10 K for the time between 
"30 and "55 min and has essentially no effect during 
other time periods (figure 7). It is interesting to 
note that although convective effects moderate the 
increase in T(r,z) max during PCM melting, total 
PCM 1 iqui f icat ion occurs earlier than for the case 
without free convection. Thus, sensible heating of 
the liquid PCM occurs and quickly increases 
T(r,z) max to about the same value that is predicted 
for the case without free convection. 

Figure 8 illustrates the fraction of total 
canister radial heat transfer (comprised of void 
conduction, void radiation, and canister wall con- 
duction) that occurs by conduction in the two canis- 
ter side walls. This "side wall fraction’’ was 
evaluated at the radial location corresponding to 
the inner surface of the outer canister wall, R 0 “. 
The side wall fraction is a quantitative measure of 
how effectively canister wails redistribute energy 
absorbed in the outer wall and also provides a qual- 
itative indication of wall temperature sensitivity 
to void location. The side wall fraction is rela- 
tively constant at "70 percent for the case with a 
void. Large perturbations in side wall fraction 
occur for brief periods of time at ~55 , "72, and 
"85 min. The three perturbations are associated 
with two step changes in the outer wall heat flux 
boundary condition and ^he point of complete PCM 
1 iqui f icat ion , respectively. For the case without a 
void, the side wall fraction is considerably lower, 
i.e., 30 to 50 percent, due to the relatively high 
thermal conductance of the PCM when compared to 
the void. The same three perturbations exist as in 
the case with the void in addition to another per- 
turbation at ~62 min. This last perturbation is 
associated with a "7 min period (from ~55 to 
"62 min) during which a thin layer of PCM freezes 
along the outer wall. Once frozen, relatively large 
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temperature gradients exist within the PCM layer 
which reduces the side wall fraction to only a cou- 
ple percent at 62'*" min. 

4. DISCUSSION 

• Results from 2D (r,z) canister analyses show 
that the largest changes in T(r,z) raax due to free 
convection and the introduction of a void are ~-10 
and ~+20 K, respectively. These changes are an 
order-of-magni tude lower than the changes predicted 
from one-dimensional analyses [ 1 ] . Therefore, this 
confirms the postulation that the effects of free 
convection and a void on canister thermal perform- 
ance are much less pronounced in a two-dimensional 
analysis than in a one-dimensional analysis because 
of the large heat transfer contribution of canister 
side walls [l] . 

4.1 On The Effects Of Free Convection 

Although relatively large Nu numbers exist dur- 
ing the PCM melting period, the length of time in 
which they occur is short. Thus over this short 
period of time, canister thermal performance in a 
1-g environment, defined in terms of maximum wall 
temperature, canister temperature gradients, and 
PCM melting rate and/or PCM utilization, is not 
greatly different from that expected in a micro-g 
environment. For canisters in the receiver cavity 
containing high quality two-phase PCM, a small 
reduction in maximum wall temperature and a small 
increase in PCM utilization would be expected 
during 1-g operation as a consequence of liquid 
PCM convection. Canisters that contain completely 
liquified PCM would experience approximately the 
same thermal cycle during 1-g and micro-g operation 
while canisters containing low quality two-phase 
PCM would be essentially unaffected. Therefore, 
within the scope of this analysis, results indicate 
that canister thermal performance during ground 
tests should not be significantly different than 
that expected on-orbit. 

The reader should be reminded that implicit in 
these results are the assumptions of (1) axisymmetry 
(which requires alignment of the gravity vector and 
the canister axis of symmetry), (2) constant radia- 
tive flux input conditions at the canister outer 
wall in each case considered, and (3) a prescribed 
void shape and location. Assumption 1 restricts the 
validity of analytical results to a ground test con- 
figuration with the canister axis vertical and with 
circumferentially uniform outer wall heating. 
Assumption 2 limits the available latitude for direct 
case-to-case comparison of results since large diff- 
erences in canister temperature predictions invali- 
dates the assumption of identical outer wall heat 
flux conditions in each case. This holds true since 
the canister outer wall boundary condition is a func- 
tion of the radiation environment in the receiver 
cavity as well as canister outer wall temperature . 
Assumption 3 removes the possibility of different 
void shapes and locations that are likely to occur 
between 1-g canister tests and micro-g canister 
operation. However, the canister performance dif- 
ferences associated with void shape and more impor- 
tantly, void location, are bounded by results from 
two cases considered in this paper, i.e., a canis- 
ter without a void and a canister with a void con- 


servatively located adjacent to the outer wall. 

This assumption is further discussed in the follow- 
ing section. 

4.2 On The Effects Of The Void 

The conservative placement of a vapor-filled 
void volume adjacent to the canister outer wall gen- 
erates the two negative effects of increased wall 
temperatures and increased wall temperature gradi- 
ents (when compared to the case without void) during 
PCM melting. These effects increase canister ther- 
mal stresses and decrease canister service life pre- 
dictions based on cumulative creep damage theory. 
Introduction of a void also reverses the sign of the 
side wall temperature gradient during the first half 
of the PCM freezing period which does not occur in 
the case without a void. The resulting change in 
the canister side wall thermal-stress distribution 
and the resulting impact (good, bad or indifferent) 
on canister life prediction is not easily determined 
without detailed structural analysis. 

A potentially beneficial effect of a void placed 
at the canister outer wall is a reduction in canis- 
ter heat loss (more precisely, canister heat 
exchange with other canisters in the receiver cavity 
and heat loss out the cavity aperture) during PCM 
freezing due to the insulating quality of the void. 
The greatest canister heat loss occurs in the hot- 
test canisters located near the aperture end of the 
receiver cavity (figure 1). These canisters are 
also located on the coolant tube near the inlet 
manifold and are thus cooled by relatively low tem- 
perature heat engine working fluid. Reduction of 
heat loss from these hottest canisters permits 
greater heat transfer to the relatively cool work- 
ing fluid near the inlet end of the coolant tube 
thereby decreasing the required working fluid heat 
transfer of canisters further down-stream. This, 
in turn, reduces the required temperature of down- 
stream canisters which must transfer heat to the 
highest temperature working fluid. 

The prescribed void shape and location were 
chosen to permit a relatively straight forward 
numerical analysis and to generate conservative pre- 
dictions of canister temperatures and stresses. To 
confirm that void heat transfer and void placement 
are conservative, consider the relationship between 
maximum canister wall temperature and the radial 
heat transfer side wall fraction. As shown in fig- 
ures 7 and 8, the canister wall radial heat transfer 
contribution can double (i.e., side wall fraction 
increases by ~2 times) with only a 20 K increase in 
maximum wall temperature during PCM melting. 

Since the side wall fraction, by definition, is 
bounded by a 100 percent value, further increases in 
canister wall heat transfer are limited to ~40 per- 
cent. Thus, regardless of the nature of void heat 
transfer, the increase in maximum canister wall tem- 
perature from introduction of a void is bounded by 
~28 K or a T(r,z) max of 1105 K for these cases. 

This T(r,z) max value is just slightly higher than 
that predicted for the canister with void case (fig- 
ure 7). Additionally, for any other viable void 
shape and location, PCM would be placed in contact 
with the outer wall thereby lowering the side wall 
fraction (figure 8, ’no void" curve) and hence, low- 
ering the maximum wall temperature. Therefore, 
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differences in void shape and location between can- 
ister 1-g ground tests and micro-g operation lead 
to relatively small changes in predicted canister 
thermal performance which are essentially bounded by 
the void cases considered in this paper. 

5. SUMMARY OF RESULTS 

Results from two-dimensional analyses of the 
PCM containment canister without a void, with a 
void, and with a void plus liquid PCM free convec- 
tion revealed several attributes of canister thermal 
performance. For engineering purposes, the change 
in canister thermal performance from 1-g to micro-g 
environments (due to free convection) was expected 
to be small if the canister ground test configura- 
tion was consistent with analytical assumptions. 

The introduction of a void and its associated poor 
heat transfer qualities increased maximum canister 
wall temperatures by ~20 K and nearly doubled canis- 
ter wall temperature gradients over those values 
predicted for the case without a void. At the same 
time, it was shown that wall temperature increases 
are bounded for any given case. It was also shown 
that for any viable void shape and location encoun- 
tered during 1-g or micro-g operation, predicted 


canister thermal performance differences were small 
and the limits of canister thermal performance were 
known . 
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FIGURE 4. - CANISTER TEMPERATURE CONTOUR (DEG K) AND PCM PHASE MAPS WITHOUT FREE CONVECTION 
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FIGURE 5. - CANISTER TEMPERATURE CONTOUR (DEG K) AND PCM PHASE MAPS WITH FREE CONVECTION. 
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FIGURE 8. - RADIAL HEAT TRANSFER SIDE WALL FRACTION. 
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